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DESCRIPTION 

SYMBOL POINT ESTIMATING APPARATUS, METHOD, PROGRAM, AND 

RECORDING MEDIUM 

TECHNICAL FIELD 

The present invention relates to estimation of symbol points of a 

signal. 

BACKGROUND ART 

There have conventionally been practiced demodulation of a received 
signal and modulation analysis of a received signal. On this occasion, it is 
necessary to precisely measure symbol points of the signal. The 
measurement of the symbol points of the received signal requires steps 
including' (l) A/D conversion of the received signal, (2) filtering to remove 
noises, (3) extraction of a symbol rate component, (4) calculation of a phase, 
and (5) conversion of the phase to a time delay (refer to a patent document 1 
(Japanese Laid Open Patent Pubhcation (Kokai) No. 2003-152816), for 
example). 

However, according to the above prior art, if the filter used to remove 
the noises has an adverse effect on firequency characteristics of the received 
signal, the symbol points of the signal cannot be precisely measured. 
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Moreover, it is necessary to carry out over sampling for the A/D conversion of 
the received signal to prevent aliasing, and there thus increase a memory 
capacity required to record the results of the A/D conversion of the received 
signal. 

It is an object of the present invention to more precisely measure 
symbol points of a received signal. 

DISCLOSURE OF THE INVENTION 

According to an aspect of the present invention, a symbol point 
estimating apparatus that estimates a symbol point of a received signal by 
determining a time delay between a sampling point of the received signal 
sampled at a sampling frequency, and the symbol point of the received signal, 
includes^ a multiphcation/sum of products output unit that outputs a sum of 
products of respective products obtained by multiplsdng a complex conjugate 
of a frequency component of an ideal signal and a frequency component of 
the received signal and a sampling angular frequency; and a time delay 
determining unit that determines a time delay to minimize an error 
component between the ideal signal and the received signal based on the 
output of the multiphcation/sum of products output unit. 

According to the thus constructed symbol point estimating apparatus, 
a symbol point estimating apparatus that estimates a symbol point of a 
received signal by determining a time delay between a sampling point of the 
received signal sampled at a sampling frequency, and the symbol point of the 
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received signal can be provided. 

The multiplication/sum of products output unit outputs a sum of 
products of respective products obtained by multiplying a complex conjugate 
of a frequency component of an ideal signal and a frequency component of 
the received signal and a sampling angular frequency. The time delay 
determining unit determines a time delay to minimize an error component 
between the ideal signal and the received signal based on the output of the 
multiplication/sum of products output unit. 

According to the present invention, the multipHcation/sum of 
products output unit may include- a frequency component product output 
unit that outputs the product of the complex conjugate of the frequency 
component of the ideal signal and the frequency component of the received 
signal; and a sum of products output unit that outputs the sum of products of 
the respective outputs of the frequency component product output unit and 
the samphng angular frequency. 

According to the present invention, the frequency component product 
output unit may include- an ideal signal frequency component output unit 
that outputs the frequency component of the ideal signal; a received signal 
frequency component output unit that outputs the frequency component of 
the received signal; a complex conjugate output unit that outputs the 
complex conjugate of the output of the ideal signal frequency component 
output imit; and a frequency component product output unit that multiplies 
the output of the complex conjugate output unit and the output of the 
received signal frequency component output unit by each other, and then 
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outputs a result of the multiplication. 

According to the present invention, the frequency component product 
output unit may include- a convolution output unit that outputs a 
convolution of the complex conjugate of the ideal signal and the received 
signal; and a frequency component output unit that outputs a frequency 
component of the output of the convolution output unit. 

According to the present invention, the sum of products output unit 
may include" a real part sum of products output unit that outputs a sum of 
products of the real part of the respective outputs of the frequency 
component product output unit and the sampling angular frequency; an 
imaginary part sum of products output unit that outputs a sum of products 
of the imaginary part of the respective outputs of the frequency component 
product output xmit and the sampling angular frequency; and a complex 
number output unit that outputs a complex number whose real part is the 
output of the real part sum of products output unit and whose imaginary 
part is the output of the imaginary part sum of products output xmit. 

According to the present invention, the time delay determining unit 
may determine the time delay based on the argument of the output of the 
multipHcation/sum of products output unit, the sampling angular frequency, 
and an error calculation length which is the number of the components of the 
received signal used to calculate the error component. 

According to the present invention, the time delay determining unit 
may include- an argument output unit that receives the output of the 
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multiplication/sum of products output unit, and outputs the argument 
thereofi and a time delay calculating unit that calciilates the time delay 
based on the output of the argument output unit, the sampling angular 
frequency, and the error calculation length. 

Another aspect of the present invention is a symbol point estimating 
method that estimates a symbol point of a received signal by determining a 
time delay between a samphng point of the received signal sampled at a 
sampling frequency, and the symbol point of the received signal, including: a 
multiplication/sum of products output step of outputting a sum of products of 
respective products obtained by multiplsdng a complex conjugate of a 
frequency component of an ideal signal and a frequency component of the 
received signal and a sampling angular frequency; and a time delay 
determining step of determining a time delay to minimize an error 
component between the ideal signal and the received signal based on the 
output of the multipHcation/sum of products output step. 

Another aspect of the present invention is a program of instructions 
for execution by the computer to perform a symbol point estimating process 
that estimates a symbol point of a received signal by determining a time 
delay between a sampling point of the received signal sampled at a samphng 
frequency, and the symbol point of the received signal, the symbol point 
estimating process including: a multiplication/sum of products output step of 
outputting a sum of products of respective products obtained by multiplying 
a complex conjugate of a frequency component of an ideal signal and a 
frequency component of the received signal and a sampling angular 
frequency; and a time delay determining step of determining a time delay to 



5 



1^ 



FADT0413PCT 



minimize an error component between the ideal signal and the received 
signal based on the output of the multiphcation/sum of products output step. 

Another aspect of the present invention is a computer-readable 
medium having a program of instructions for execution by the computer to 
perform a symbol point estimating process that estimates a symbol point of a 
received signal by determining a time delay between a sampHng point of the 
received signal sampled at a sampling frequency, and the symbol point of the 
received signal, the symbol point estimating process including- a 
multiplication/sum of products output step of outputting a sum of products of 
respective products obtained by multiplying a complex conjugate of a 
frequency component of an ideal signal and a frequency component of the 
received signal and a sampling angular frequency; and a time delay 
determining step of determining a time delay to minimize an error 
component between the ideal signal and the received signal based on the 
output of the multiplication/sum of products output step. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a configuration of a symbol point 
estimating apparatus 1 according an embodiment of the present invention 

FIG. 2 is a diagram showing an EVM, which is an error component 
between an ideal signal r(k) and a received signal z(k); 

FIG. 3 is a diagram showing a configuration of a frequency 
component product output imit 121 

FIG. 4 is a diagram showing a variation of the configuration of the 
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frequency component product output unit 12; 

FIG. 5 is a diagram showing a configuration of a sum of products 
output unit 14; and 

FIG. 6 is a diagram showing a configuration of a time delay 
determining unit 20. 

BEST MODE FOR CARRYING OUT THE INVENTION 

A description will now be given of a best mode to carry out the 
present invention with reference to drawings. 

FIG. 1 is a block diagram showing a configuration of a symbol point 
estimating apparatus 1 according an embodiment of the present invention. 
The symbol point estimating apparatus 1 is used to estimate symbol points of 
a received signal z(k). The estimation of the symbol points enables 
demodulation of the received signal z(k) and modulation analysis of the 
received signal z(k). The estimation of the symbol points of the received 
signal z(k) is carried out by determining a time delay T between sampling 
points of the received signal z(k) sampled at a sampHng frequency fs and the 
symbol points of the received signal z(k). 

The symbol point estimating apparatus 1 includes a 
multipHcation/sum of products output unit 10 and a time delay determining 
unit 20. 

The mvdtipHcation/sum of products output unit 10 outputs a sum of 
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products Aei ^ of respective products Y(n)=Z(n)R(n)* obtained by multiplying 
a complex conjugate R(n)* of a frequency component R(n) of an ideal signal 
r(k) and a frequency component Z(n) of the received signal z(k) by each other, 
and a sampling angular frequency A O) (=2 TVfsfN). It should be noted that 
N denotes an EVM calculation length. Moreover, the ideal signal r(k) is 
generated from the received signal z(k). It should be noted that EVM (Error 
Vector Magnitude) is an error component between the ideal signal r(k) and 
the received signal zCk) as shown in FIG. 2. The EVM is defined by the 
following equation (l). It should be noted that N denotes the EVM 
calculation length. 



The multiplication/sum of products output unit 10 includes a 
fi-equency component product output unit 12, and a sum of products output 
unit 14. 

The firequency component product output unit 12 outputs a product 
Y(n)=Z(n)R(n)* of a complex conjugate R(n)* of a frequency component R(n) 
of the ideal signal r(k) and the frequency component Z(n) of the received 
signal z(k). A configuration of the frequency component product output unit 
12 is shown in FIG. 3. The frequency component product output unit 12 
includes an FFT unit Gdeal signal frequency component output means) 122, 



[EQU. 1] 




[%rms] 



(1) 
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an FFT unit (received signal frequency component output means) 124, a 
complex conjugate output unit 126, and a multiplier (frequency component 
product output means) 128. 

The FFT unit Gdeal signal frequency component output means) 122 
applies the FFT (Fast Fourier Transform) to the ideal signal. r(k), and 
outputs a result thereof. The result of the FFT applied to the ideal signal 
r(k) is the frequency component R(n) of the ideal signal r(k). 

The FFT unit (received signal frequency component output means) 
124 appHes the FFT (Fast Fourier Transform) to the received signal z(k), and 
outputs a result thereof. The result of the FFT applied to the received 
signal z(k) is the frequency component Z(n) of the received signal z(k). 

The complex conjugate output unit 126 outputs the complex 
conjugate R(n)* of the output R(n) of the FFT unit Gdeal signal frequency 
component output means) 122. 

The midtiplier (frequency component product output means) 128 
multiplies the output R(n)* of the complex conjugate output unit 126 and the 
output Z(n) of the FFT unit (received signal frequency component output 
means) 124 by each other, and outputs a result thereof. This output is 
Y(n)=Z(n)R(n)*. 

A variation of the configuration of the frequency component product 
output unit 12 is shown in FIG. 4. As shown in FIG. 4, the frequency 
component product output unit 12 includes a complex conjugate output unit 
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121, a convolution output unit 123, and an FFT unit (frequency component 
output means) 125. 

The complex conjugate output unit 121 outputs the complex 
conjugate r(k)* of the ideal signal r(k). 

The convolution output unit 123 outputs a convolution of the output 
r(k)* of the complex conjugate output unit 12 1 and the received signal z(k). 

The FFT unit (frequency component output means) 125 applies the 
FFT (Fast Fourier Transform) to the output of the convolution output unit 
123, and outputs a result thereof. The result of appljdng the FFT to the 
output of the convolution output imit 123 is Y(n)=Z(n)R(n)*. 

The sum of products output unit 14 outputs a sum of products Ae* ^ of 
the output Y(n) of the frequency component product output unit 12 and the 
sampUng angular frequency A O). 

A configuration of the sum of products output unit 14 is shown in FIG. 
5. The sum of products output unit 14 includes a real part acquisition unit 
141, a real part sum of products calculation unit 142, an imaginary part 
acquisition unit 143, an imaginary part sum of products calculation unit 144, 
and a complex number output unit 146. 

The real part acquisition unit 141 acquires the real part I(n) of Y(n). 

The real part sum of products calculation unit 142 outputs a sum of 
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products of I(n) and the sampling angular frequency A oj. The sum of the 
products of I(n) and Ao) is represented by the following equation (2). 

[EQU. 2] 

AT/ 2-1 

y|n|A<y/(«) (2) 

The imaginary part acquisition unit 143 acquires the imaginary part 
Q(n) of Y(n). 

The imaginary part sum of products calculation unit 144 outputs a 
sum of products of Q(n) and the samphng angular frequency A w . The sum 
of the products of Q(n) and A 6c> is represented by the following equation (3). 

[EQU. 3] 

AT/ 2-1 

J«Ao;e(«) (3) 

The complex number output unit 146 outputs a complex number 
whose real part is the output of the real part sum of products calculation imit 
142 and whose imaginary part is the output of the imaginary part sum of 
products calculation unit 144. The output of the complex number output 
unit 146 is represented as Ae^^. The complex number output unit 146 
includes a multiplier 146a and an adder 146b. The multipHer 146a 
multipHes the output of the imaginary part sum of products calculation unit 
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144 by j to obtain an imaginary number. The adder 146b adds an 

output of the multipUer 146a to the output of the real part sum of products 
calculation unit 142. The output of the adder 146b is Ae*^. 

The time delay determining unit 20 determines the time delay T so 
as to minimize the error component (EVM) between the ideal signal r(k) and 
the received signal z(k) based on the output Aei^ of the multiphcation/sum of 
products output unit 10. 

The EVM is obtained by normalizing and then extracting the square 
root of an error component s defined by the following equation (4). 

[EQU. 4] 

Thus, the error component a is minimized to minimize the EVM. 
When the error component £ is minimized, the following equation (5) holds. 
Namely, a partial derivative of the error component a with respect to the 
time delay T is 0 (zero). 



N-l 



z{k-ry{k) 



(4) 
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[EQU.5] 



— =0 (5) 



There thus can be obtained a time delay T which minimize the 
EVM by solving the equation (5) for the time delay T. An equation (6) is 
obtained by solving the equation (5) for the time delay T (proof is provided 
later). 

[EQU. 6] 
40 

^'AwN (6) 

Thus, the time delay determining unit 20 can determine the time 
delay T based on the argument 0 of the output A& ^ of the 
multipUcation/sum of products output unit 10, the sampling angular 
frequency A O), and the EVM calculation length N. 

A configuration of the time delay determining unit 20 is shown in 
FIG. 6. The time delay determining unit 20 includes an argument output 
unit 22 and a time delay calculating unit 24. 

The argument output unit 22 receives the output A& ^ of the 
m\dtiplication/sum of products output unit 10, and outputs the argument 0 
thereof. The time delay calculating unit 24 calculates the time delay T 
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based on the output Q of the argument output unit 22, the sampling 
angular frequency A CO, and the error calculation length N. Specifically, 
the time delay T is calculated by assigning 0 , b^uy, and N to the right 
side of the equation (6). The time delay T determined in this way makes 
the equation (5) hold, and thus minimizes the error component e . Thus, 
the error component (EVM) can be minimized. : 

A description will now be given of an operation of the embodiment of 
the present invention. 

First, the ideal signal r(k) is generated from the received signal z(k). 
The received signal z(k) and the ideal signal r(k) are supplied to the 
frequency component product output unit 12 of the multiplication/sum of 
products output unit 10. The frequency component product output unit 12 
outputs Y(n)=Z(n)R(n)*. The sum of products output unit 14 obtains the 
sum of the products of Y(n) and the sampUng angular frequency A w, and 
outputs the result as Aei ^ . 

The resulting sum of products Ae^^ is supplied to the time delay 
determining unit 20. The time delay determining unit 20 calculates the 
time delay T based on the argument Q of Aei^, the sampling angular 
frequency Act), and the EVM calculation length N. The determined time 
delay T can minimize EVM. 

According to the present embodiment, it is possible to determine the 
time delay T according to the frequency components (Z(n), R(n)) of the 
received signal z(k) and the ideal signal r(k). Then, it is possible to estimate 
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the symbol points of the received signal z(k) according to the time delay T . 
Since the frequency components (Z(n), R(n)) are used on this occasion, it is 
possible to more precisely estimate the symbol points of the received signal 
z(k) compared with the conventional case where the temporal components 
(z(k), r(k)) are used. 

Moreover, the above -described embodiment may be realized in the 
following manner. Namely, a computer is provided with a CPU, a hard disk, 
and a media (such as a floppy disk (registered trade mark) and a CD-ROM) 
reader, and the media reader is caused to read a medium recording a 
program realizing the above-described respective parts (such as the 
multiplication/sum of products output unit 10 and the time delay 
determining unit 20), thereby installing the program on the hard disk. This 
method may also realize the above- described functions. 

[Proof of Obtainment of Equation (6) from Equation (5)] 

First, the error component a is represented by the frequency 
component R(n) of the ideal signal, and the frequency component Z(n) of the 
received signal. 

First, discrete Fourier transform pairs of z and r are represented as- 
z(k)<^Z(n) and r(k)<^R(n). On this occasion, the following equation (7) holds 
according to Parseval's equality. 
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[EQU. 7] 



> (7) 



K^K(A:) = ^|z(«K(«) 



On this occasion, for Z(n) and RGO, Z(n)=ZGi-N) and R(n)=R(n-N) 
hold, and the equation (7) is thus rewritten as the following equation (8). 



[EQU. 8] 



N-l ^ 1 N/2-1 



N-l ^ 1 N/2-1 



V (8) 



N-l 1 N/2-1 

2z{ky{k) = j- fz{n)R'{n) 



J 
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Moreover the following equation (9) holds according to the time 
shifting theorem. 

[EQU. 9] 

z{k - r) e-'"'^'^z{n) ( g ) 

When the equations (8) and (9) are assigned to the equation (4) which 
defines the error component e , the following equation (10) is obtained. 

[EQU, 10] 



N-l 



e = y\z{k -xf + y |r(A:)|' - 2Re 



N/2-1 



r N/2-i 



jr "fUnf I'm' -|Re Ve-'--'Z(„)«-(») 



r N/2-i 



Nf2^\ A?72-l 



N/2-l 



(10) 



Then, a third term of the equation (lO) is transformed to obtain the 
following equation (ll). 
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[EQU. 11] 



r Nii-i 



Re 



AT/ 2-1 



(1 1) 



The equation (ll) represents a relationship among terms of the real 
part, and does not include terms of the imaginary part. Thus, even for 
complex conjugates of respective terms of the equation (ll), the equation (ll) 
still holds. The following equation (12) is thus obtained by replacing the 
first term on the right side of the equation (ll) by the complex conjugate 
thereof, and then transforming the equation (ll). 

[EQU. 12] 



Re 



■ N/2-1 


= Re 







-1 f X* N/2-1 



= Re 



N/2-l 



TNf2 



N/2-1 



= Re Y e-^"^^y * (- Ai) + y e-'^'^Yin) 



(12) 



When the equation (12) is assigned to the equation (lO), the error 
component a is represented by the following equation (l3): 



18 



FADT0413PCT 



[EQU. 13] 



1 N/2-1 ^ 1 iV/2-1 ^ ^ VN/1 N/2-i 

(13) 



There is then obtained the value (left side of the equation (5)) which 
is the partial derivative of the error component e with respect to the time 
delay T. 

The value obtained by partially differentiating the error component 
£ by the time delay T is represented by the following equation (14). 

[EQU. 14] 



a£ 

dr ~ N 



a Re 



JV/2 JV/2-1 



dr 



7 rAr/2 

= -— Re ye-"^"" 



N/2-l 



(- jnAcoy (- n) + ^ e-^^"^ (- ;« A«>)k(/j) 



9 rN/2 N/2-1 



(14) 



On this occasion, the time delay T is small to a certain extent, and a 
relationship represented by the following equation (15) holds for n=0 to N/2. 
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[EQU. 15] 



4 



(15) 



The following equation (16) is obtained by assigning the equation (14) 
to the equation (15). 

[EQU. 16] 



ds __2_ 
dr ~~N 



Im 



A I 



(N/2 N/2-1 

J Y «Aa>y *(-«)+ ynAcoY{n) 



Im 

N 



j—Acm r N/2 Y 

* \ j|n|Acy/(n)+7 j«Aa>i2(«) 



(16) 



According to the definition of the sum of products output unit 14, the 
following equation (17) holds. 

[EQU. 17] 



N/2 N/2 

y\n\A(ol{n)+ j VnAcoQin) ^ Ae'^ (17) 



The equation (16) can thus be rewritten as the following equation 

(18). 
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[EQU. 18] 



dr N 



N 




(18) 



There is finally to be obtained such T that a partisil derivative of 
the error component £ with respect to the time delay t is 0. 

It is necessary that the following equation (19) holds in order to cause 
the equation (18) to be 0. This is because the imaginary part is 0 if the 
argument of a complex number is 0°. 



The following equation (20) is obtained by solving the equation (19) 
with respect to T . 



[EQU. 19] 



—Aajt -d = 0 
4 



(19) 



[EQU. 20] 



T = 



AwN 



(2 0) 
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The equation (20) is the same as the equation (6). 

Thus, it is possible to obtain the equation (6) from the equation (5). 
[End of Proof] 
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